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ABSTRACT 

We present the results of a search for bright high-redshift galaxies in two large SCUBA 
scan-maps of Galactic regions. A Mexican Hat Wavelet technique was used to locate 
point sources in these maps, which suffer high foreground contamination as well as 
typical scan-map noise signatures. A catalogue of point source objects was selected and 
observed again in the submillimetre continuum, and in HCO + (3— >2) at zero redshift 
to rule out Galactic sources. No extragalactic sources were found. Simulations show 
that the survey was sensitive to sources with fluxes > 50 mJy, depending on the local 
background. These simulations result in upper limits on the 850-/mi counts of SCUBA 
galaxies of 53 per square degree at 50 mJy and 2.9 per square degree at 100 mJy. 
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1 INTRODUCTION 

The Submill imetre Common -User Bolometer Array 
(SCUBA, IHolland et alJ Il999l) has been operating at 
the James Clerk Maxwell Telescope (JCMT) since 1997, 
and in that time many blank-field extragalactic sur- 
veys have been c ompleted, as summarised in Table [T 
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SCUBA galaxies over the flux range 2-15 mJy are now well 
measured. This is illustrated in Fig. where the counts 
from several surveys have been compiled. 

This range of fluxes is imposed both by the abilities of 
SCUBA and the JCMT, and by the nature of the counts. 
SCUBA has a beam at 850 /jm of about 14", which, in com- 
bination with the expected slope of the counts and the clus- 
tering strength of these objects, results in a confusion limit 
for reliable detections of ab out 2 mJy iBlain et alJll998l: 
lHogdl200ll : ISmail et al.ll2002l) . 

At the higher flux end, the limit of ~ 15 mJy is set 
by the steep slope of the counts, and by the sensitivity of 
SCUBA. The variation of the integral galaxy count per unit 
area, N, with limiting flux density, S, can be described as 




Flux (mJy) 

Figure 1. The integral counts of SCUBA galaxies as measured 
by the surveys listed in Table HI The solid line shows a power law 
fitted to the bright count end so that N oc S~ 2 where N is the 
number of galaxies per sky and S is the 850/im flux. 



a power law, N oc S~ a . Between 5 and 15 mJy, a has 
been constrained at ~ 2 (Fig. 0. Beyond about 15 mJy 
the counts are not well known, although at some brighter 
flux they are expected to tend to the low-redshift Euclidean 
slope with a — 1.5. The recent detection o f a bright sub- 
millimetre quasar, with Ss50ftm ~ 25 mJy (iKnudsen et alJ 
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Table 1. A summary of the blank field surveys carried out using SCUBA. Lensing survey sensitivities have been corrected for lensing. 
Symbol in the fifth column refers to that used in Fig. 



Survey 


Depth (3ct, mjy) 


Area (square arcmin) 


Reference 


Symbol 


Hawaiian deep fields 
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Bareer et al. 1999 


Filled triangle 
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£C 






8 


104 


tl 




HDF jiggle-map 


1.5 
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Huehes et al. 1998 


Filled star 
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Eales et al. 2000 


Crossed circle 
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Webb et al. 2002 
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Borvs et al. 2002 
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Canada Lens Survey 
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2003), strengthens the idea that there is some bright, as yet 
undiscovered, population of SCUBA galaxies. 

In this paper, we address the question of the very bright 
count by searching for serendipitous observations of external 
galaxies in observations of Galactic fields. Two of the largest 
SCUBA maps made so far ar e of the Galactic centre a nd the 
Perseus star-forming regions l)Pierce-Price et, al.l200fi Qual- 
trough et al., in prep.), which total approximately 5 square 
degrees. SCUBA scan-maps pointed towards known Galac- 
tic features are the only SCUBA maps on such large scales. 
Whilst each map contains bright Galactic emission, this is 
of course optically thin in the submillimetre and hence any 
background extragalactic objects will shine through. Hence 
the challenge is to separate the two astronomical compo- 
nents. Furthermore, SCUBA scan-maps are prone to com- 
plicated noise structures, which arise both from the usual 
atmosphere, telescope and instrument sources and from the 
image reconstruction technique. Thus a source extraction 
routine which does not rely on assumptions about the char- 
acteristics of the noise (astronomical and otherwise) in a 
map is preferable. For this purpose, a Mexican Hat Wavelet 
(MHW) algorithm was used. 



2 SCAN-MAPPING WITH SCUBA: THE DATA 

The scan-mapping mode available for SCUBA observations 
is the principal way to make large, shallow maps with 
SCUBA. Relatively rarely used in comparison to the jiggle- 
map technique, scan-maps thus far have been predominantly 
of Galactic regions. The only completely blan k-field scan- 
map i s the survey of the Hubble Deep Field bv lBorvs et alJ 
J2002D . 

In scan-map mode <|,Tenness et al.ll2000l L the SCUBA 
bolometer array is scanned across the sky to cover an area 
of approximately 10' x 10', much larger than the single- 
pointing 2.3'diameter field of view. Whilst scanning, a 'chop' 
movement at a rate of ~ 8 Hz is used to provide a reference 
measurement of atmospheric emission. The secondary mir- 
ror of the JCMT switches the array between two patches of 
sky separated by the chop throw, which is sufficiently small 
to enable accurate removal of near-field sky emission for ob- 
jects smaller than the chop throw scale, such as the point 
sources we are interested in. 



SCUBA scans each field six times, using three chop 
throws each in two chop directions, chosen so that iden- 
tical chop throws are at 90° to each other. The six scans 
are then added together by sensitivity-weighted averaging 
within Fourier space to produce the final image. However, 
the coverage of the Fourier frequency plane is imperfect, re- 
sulting most noticeably in low sensitivity of the array to 
the larger spatial scales within the overall scan field. The 
14"beam size of the telescope also acts to reduce sensitivity 
at very small spatial scales. 

Prior to the full combination each scan goes through 
several reduction steps, common to most SCUBA data re- 
duction. The maps are flatfielded to correct for different 
sensitivities between bolometers, the extinction of the sig- 
nal due to the sky is corrected for, spikes caused by cosmic 
rays are removed, and remaining instrumental baselines and 
particularly noisy bolometers are removed. Finally, the time 
series data is rebinned onto an image grid, here using a pixel 
scale of 3"at 850 /im, and the rebinned maps are calibrated 
using maps of planets or other well-known objects taken in 
the same observing mode. These steps were all performed 
using the SCUBA User Re duction Facility software (SURF, 
iJenness fc LightfoodlT99sl) . 

In the rest of this paper, the term 'background' will 
be used to describe every contribution to the observed flux 
in images except the point sources i.e. the 'background' in- 
cludes the diffuse Galactic emission, the atmosphere, instru- 
ment and telescope noise, and any artifacts such as noise 
spikes or features related to the reconstruction of the maps 
in Fourier space. 



3 POINT SOURCE EXTRACTION 
3.1 SExtractor 

Our initial attempts to locate point sources in these maps 
used the S tarlink source extraction package known as SEx- 
tractor feertin fc Arnoutsl [l996). It is intended primarily 
for optical images, with perhaps thousands of objects which 
may be slightly or completely resolved. SExtractor thus 
identifies sources as clusters of pixels all with flux above 
some threshold, which is chosen with respect to an estimated 
background. Filters can be used to guide the required shape 
of the clustered pixels. 
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Initial runs with SExtractor found that in the 
crowded Galactic centre fields it was an unreliable detector 
of even the most obvious point sources: any source detection 
software which relies on being able to accurately characterise 
the background will not perform well with the complicated 
backgrounds found in SCUBA scan-maps. As explained be- 
low, wavelet routines do not need prior assumptions about 
the background emission in a map in order to find point 
sources, and so this detection technique was chosen. A fur- 
ther advantage of the MHW technique is that it can be finely 
tuned to finding sources of the appropriate size, leading to 
fewer false positives. 




3.2 The Mexican Hat Wavelet 

The MHW technique has become a commonly-used tech- 
nique to separate the point sources from future Planck data 
JCavon et alJ l2OO0t ISanz et al.l boOlt IVielva et alj l2001bl 
2003). The combination of the MHW techniq ue with the 
Maximum Entropy method jHobson et alJl99Sf) has proved 
to be very successful in se parating all the com ponents of the 
simulated microwave sky (IVielva et alj|2001al) . 

Other wavelet techniques have been sug gested for ex- 
tract ion o f point sources i n CMB maps llTenorio et alJ 
I1999T) . but ISanz et alJ feOOlf) demonstrated that for Gaus- 
sian sources the best wavelet basis to use for a wide range 
of backgrounds is the Mexican Hat. This was also found 
to be preferable to a straightforward Gaussian filter for all 
backgrounds except the (unrealistic) case of white noise. A 
Gaussian linear filter generally produces a smoothing effect, 
lowering the amplitudes of detected sources with respect to 
the noise levels within a map. 

A wavelet filter acts as a band-pass filter on the Fourier 
transform (FT) of an image - depending on the definition of 
the 'mother wavelet' (see section l4"TI . the power from certain 
spatial frequencies is amplified. In CMB maps, the removal 
of several types of foreground signal, including extragalactic 
point sources, is essential to separate out the cosmological 
signal. In the application of the wavelet algorithms in this 
paper, the extragalactic objects become the final target. The 
MHW tech nique was also used to detect point sources in X- 
ray data bv lDamiani et all il997f) . We have also applied the 
MHW technique to a jiggle-map survey of cluster lenses. The 
details of the method, which differs slightly from that used 
here for scan-maps, will be described in K. K. Knudsen et 
al., in prep. 

In Fig. [5] a scan-map image of the point source calibra- 
tor Uranus at 850 /im is shown. The contours demonstrate 
that the Gaussian peak of the beamshape is affected by ar- 
tifacts caused by the incomplete and variable coverage of 
scanning of the Fourier plane as well as imperfections in the 
surface of the JCMT's dish. These artifacts are known to be 
more prevalent in 450-/im SCUBA maps, and so the search 
for point source candidates was carried out using the 850- 
fixa maps only. Furthermore, the ratio of Galactic flux to 
typical (anticipated) extragalactic point source flux is lower 
at 850 ^m, making this the more sensitive wavelength. 



Figure 2. Scan- map images of the point source calibrator Uranus 
at 850 fim as observed in f999, coeval with the scan-map images 
of the Galactic centre used in this paper. The contours begin at 
4 Jy beam -1 . 



4 MATHEMATICS AND APPLICATION OF 
THE MEXICAN HAT WAVELET 

4.1 Mathematics 

In this section the standard mathematics involved in the 
Mexican Hat transform will be described. Further detail 
can be found in the original papers d escribing this work 
JCavon et alJl200d IVielva et al.ll2001lJl . 

The continuous isotropic wavelet transform of a 2- 
dimensional signal f(x) is 



w f (R,b) = 



\x — b\ 
R 



/(x), 



(1) 



where u>f(R, b) is the wavelet coefficient associated with the 
scale R at the point b. ip jj-^ is known as the 'mother 
wavelet' which defines the spatial frequencies of interest. The 
two-dimensional Mexican Hat mother wavelet is given by 



2 /2R 2 



(2) 



A Gaussian point source, generated by the convolution 
of a delta function of amplitude A with a telescope beam of 
FWHM 8, can be described as 



-x 2 /2B 2 



(3) 



where Q is the area under the beam. For such an input 
signal, the wavelet coefficient for a scale R at the position 
of the source is 



R 



= 2V^ 



[R/ef 



n [i + {R/Q) 2 



(4) 



The effectiveness of the Mexican Hat Wavelet is re- 
vealed through the definition of the detectability or detection 
level of the source. In real space, the criterion for detection 
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of a source is usually the signal-to-noise of the source, given 
by 



D r = 



A/n 



(•») 



where cr n is the dispersion of the real space noise. 

In wavelet space, the equivalent parameter for the same 
point source (as a function of scale R) is defined as 



D„(R) 



MRS 



(6) 



where a Wn (R) is the dispersion of the noise field in wavelet 
space. 

Looking at the ratio of these two detectabilities, the 
amplification, A, due to the transform into wavelet space is 



A : 



D W (R) 

D r 



= 2V2^R 



(R/ef 



[1 + (i?/6>) 2 ] 2 a Wn {R)' 



(7) 



Thus, for a suitable scale -R op t — 9, the amplification 
possible by transforming an image into wavelet space can 
be very large, because at this scale the dispersion of the 
general background in wavelet space will be lower than in 
real space. Unless the background has Gaussian peaks of the 
same scale as the point sources (i.e. unless the background 
is completely indistinguishable from the real point sources) , 
it will always be relatively less amplified. 



4.2 Candidate selection process 

As can be seen from Fig.H the scan-map technique produces 
very irregularly shaped maps. However, the MHW routines 
require square input images with no blank pixels. These were 
created by trimming each 10' x 10'map individually to keep 
the largest possible areas available for searching for point 
sources. In fact the original maps overlapped in order to 
mosaic them into the final large maps, and so the overall 
area lost by trimming the original irregular maps was very 
small. Two typical squares from this stage of the process are 
shown in Fig. [3] 

Simulations performed using the Mexican Hat w avelet 
algorithm for use wit h Planck dCavon et alJ l2000t 
IVielva et al]l2001bl I2003F) were used to determine the de- 
tection parameters for this study, since the effects of vari- 
ous background components at an identical wavelength were 
considered. Detection of a source is controlled by two param- 
eters. The algorithm for source selection was as follows: 

(i) Firstly the optimal scale, R op t, was calculated by the 
MHW software for each input map. This involved iterating 
through small changes in the value of R around the point 
source scale 6 (the beamsize) until the maximum amplifi- 
cation A was found for the map. The final value of _R op t 
depen ds upon the impact of the background on varying 
scales JVielva et al.ll2001bft - a background with a charac- 
teristic scale a little larger than 8, for instance, can be most 
strongly counteracted by a value of R op t < 0. As expected, 
in the Perseus maps, where the most prominent background 
component is on small scales (Fig. |3J, the value of _R op t 
was typically larger than beam size, with a modal value of 
1.2 x beam size. Conversely, in the Galactic centre maps, the 
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Figure 3. Two squares used as input to MHW point source detec- 
tion. The top square is a Perseus field, whilst the bottom square 
is a Galactic centre field. These two squares demonstrate the va- 
riety of backgrounds that the MHW algorithms had to cope with. 
The small-scale (noise) structure visible in the Perseus image is 
also present in the Galactic centre data, but is masked by the 
much stronger Galactic emission. 



strong large-scale Galactic emission led to a modal value of 
R op t = 0.7 x beam size. 

(ii) Point source candidates were then selected at posi- 
tions with wavelet coefficient values Wf(R opt ) > 5a Wrl (R op t) 
i.e. D w (R opt ) > 5. 

(iii) For each candidate, the 'experimental' Wf(R) was 
compared to the theoretical variation expected with R, as a 
further check on the source's shape. A value of \ was cal- 
culated between the expected and experimental results and 
a source was retained as a candidate if \ 2 < 4. 

The 339 candidates (mostly in the Galactic centre 
fields) flagged by the MHW package were each examined 
by eye, to rule out those sources which were obviously asso- 
ciated with strong filamentary Galactic features. This ruled 
out the majority of the detected point sources. Next, the de- 
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tected position of each source was tested against the SIMBAD 
online database 1 to rule out known Galactic point sources, 
as well as to check whether any known extragalactic objects 
could be identified with the submillimetre sources. None of 
the latter were located. Candidates were ruled out from fur- 
ther investigation whenever they were close, with respect 
to the angular resolution both of SCUBA and other instru- 
ments, to known Galactic objects. Typical Galactic objects 
which could appear point-like to the SCUBA 14" 850-/xm 
beam include Herbig-Haro objects or protostars, as well as 
knots of dust within larger structures, especially a long fila- 
ments JPierce-Price et aljl200ot IVisser et alJl2002l) . 

From this process a final list of 21 candidates was ob- 
tained, 15 in the Perseus region and 6 in the Galactic centre 
region 2 . Of the original list, 94 per cent were rejected either 
by association with clear Galactic features or by the SIMBAD 
comparison. The 21 candidates, with their MHW-estimated 
fluxes at 850 (jm, are listed in table 



5 FOLLOW-UP OBSERVATIONS OF THE 
CANDIDATE SOURCES 

To verify the point source candidates found by MHW, sec- 
ondary observations of their locations with SCUBA were 
carried out in August 2001. The aim of these observations 
was to: i) verify their point source nature, by re-imaging 
them in jiggle-map mode at both 850 and 450 /im - the 
smaller beam of the 450-//m observations providing a check 
that the candidates, if detected at 450 /im, still appeared 
point-like on smaller scales; ii) determine accurate fluxes and 
positions for the candidates, in a map with a single pointing 
at each location. 

5.1 SCUBA jiggle-maps of the point source 
candidates 

During the nights of the 26th, 30th and 31st of August 
2001, each of the 21 candidate positions was observed with 
SCUBA. Each position was observed for long enough that, 
based on the estimated fluxes from the MHW output, the 
candidates should be detected at more than 5-cr levels in 
both wavebands. The 450- /im fluxes were estimated assum- 
ing a spectral index in the submillimetre of f3 ~ 1.5, where 
flux S u (xv 13 . Details of these observations are listed in table 

m 

Standard SCUBA jiggle-maps produce a single 2.3' 
field, centred on the candidate position. This mode is so- 
called because during observations the array is 'jiggled' to 
ensure the entire image is fully sampled. This is unnecessary 

1 http:/ /simbad.u-strasbg.fr/Simbad 

2 Included in the Galactic centre list was an object (gc6) which, 
whilst only weakly detected by the MHW routine, had extremely 
strong fluxes and was apparently point-like at both 850 and 450 
/im. This object was retained because despite its bright flux, no 
other survey in SIMBAD showed a detection nearby. Its bright flux 
indicated that only a very short integration time would be neces- 
sary to confirm it as a real source. The reason for its very weak 
detection by the wavelets routines despite its shape and strength 
was not clear at this point (but became clear following later sim- 
ulations - see section [6.3l . 



in scan-map mode as the array moves across the sky and so 
fills in the gaps between bolometers. In jiggle-map mode, as 
well as chopping, an extra 'nod' movement is included. After 
chopping between the on-source position and an off-source 
position to one side of the source, the telescope nods so that 
the secondary mirror now chops between the on-source posi- 
tion and a new off-source position, at the same distance but 
on the opposite side of the on-source position. The main 
advantage of using the nod and chop combination is that 
an estimate of any overall gradient in the sky emission can 
be made. The chop throws and directions were chosen indi- 
vidually for each source to avoid chopping onto particularly 
strong areas of Galactic emission visible in the scan-maps. 

All of the maps from the jiggle-map observations 
were reduced using stan dard Surf reduction routines 
jjenness fc Lightfoodll998l) . following the same steps as de- 
scribed above for scan-maps: flatfielding, extinction correc- 
tion, spike removal, noisy bolometer removal, and rebinning. 
Finally all were calibrated using planets or other well-known 
point sources. 

The success rate of locating the point source candidates 
was low. Three maps, corresponding to the sources labelled 
gel, gc2 and gc6 in Table [5] do show point source objects 
at their centre, whereas the remainder show either much 
larger-scale structure or nothing at all. The maps of these 
three sources are shown in Fig. [I] along with an example of 
an empty field where the candidate clearly is absent. The 
450-pim data for the three detected objects are shown in 
Fig. |K| The brightest source at 850 /im, gc6, can be seen in 
the higher frequency map, but the other two cannot. None 
of the other candidates had 450-pm detections of a point- 
source nature. The details of the detections of these three 
objects in the SCUBA maps are given in table |3] 

5.2 Molecular line observations 

To verify whether the three detected point source objects 
could be extragalactic, they were observed using the het- 
erodyne Receiver A3 on the JCMT, at the Galactic centre 
rest frequency of the line HCO + (3^2). This line is a useful 
tracer of very dense molecular gas - the critical H2 density 
to populate the H CO + (J = 3) state is ~ 7.5 x 10 7 cm~ 3 
JWild et alJll992h . within the range expected for molecu- 
lar clouds responsible for star formation. HCO + (3— >2) was 
used in preference to CO lines as it is less widely distributed 
in the Galaxy, and so the likelihood of confusion with an- 
other Galactic source is very low. The absence of a Galactic 
HCO + (3— >2) line in these objects would suggest that an 
extragalactic origin for the sources could be plausible. 

The observations were carried out in frequency-switched 
mode, which is analogous to a nod and chop in frequency 
rather than in space. The signature of a detection in fre- 
quency switched mode is therefore a negative-positive- 
negative shape in frequency space. Raster maps were made 
which pointed at a square of 5 x 5 positions and covered 
an area of 40" x 40"around the MHW-derived positions, to 
allow for possible positional errors in the predictions. The 
resolution of the JCMT at the frequency of HCO + (3-^2) is 
about 18", slightly larger than the SCUBA 850- /xm beam 
size. 

The data were r educed using the stand ard routines from 
the specx package iPrestaee et aljEoobT) . The results are 
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Table 2. Observations of point source candidates from the Galactic centre and Perseus region maps. Positions (J2000) and fluxes are 
those predicted by the MHW routines. The final three columns give observing details: observation date, time spent integrating on each 
source and average sky opacity during observation. 



UU LI I Ct^, 


R. A. 


Dec. 




\JUB. Lid/lit^ 


T- * Cs'l 


r 850^m 


gel 


17 46 42.7 


-28 08 20 


319 


26/08/01 


640 


0.340 


gc2 


17 44 19.2 


-29 33 35 


344 


26/08/01 


640 


0.340 


gc3 


17 44 05.1 


-29 41 53 


224 


26/08/01 


1024 


0.340 


gc4 


17 44 13.9 


-29 40 47 


245 


26/08/01 


768 


0.340 


gc5 


17 44 21.4 


-29 49 04 


232 


26/08/01 


1024 


0.340 


gc6 


17 43 10.3 


-29 51 44 


1556 


26/08/01 


384 


0.340 


perl 


03 25 41.9 


+30 40 57 


83 


31/08/01 


1280 


0.265 


per2 


03 25 35.6 


+30 31 27 


100 


30/08/01 


1920 


0.350 


per3 


03 25 34.6 


+30 30 36 


63 


31/08/01 


2560 


0.275 


per4 


03 25 56.5 


+30 39 52 


143 


30/08/01 


1280 


0.336 


per5 


03 26 24.0 


+30 20 28 


109 


31/08/01 


1280 


0.275 


per6 


03 26 51.8 


+30 28 09 


98 


31/08/01 


1280 


0.330 


per7 


03 45 16.9 


+32 04 47 


137 


30/08/01 


1280 


0.350 


pcr8 


03 28 13.5 


+31 36 24 


158 


30/08/01 


1024 


0.336 


per9 


03 34 54.6 


+31 12 15 


84 


31/08/01 


1280 


0.330 


perlO 


03 33 04.2 


+30 51 45 


101 


31/08/01 


1280 


0.275 


peril 


03 32 40.3 


+31 03 03 


121 


31/08/01 


1280 


0.265 


perl2 


03 33 01.7 


+30 52 24 


141 


30/08/01 


1280 


0.358 


perl3 


03 33 34.1 


+31 20 09 


136 


30/08/01 


1280 


0.358 


perl4 


03 36 42.7 


+31 14 40 


86 


31/08/01 


1280 


0.330 


perl5 


03 36 40.1 


+31 05 10 


96 


31/08/01 


1280 


0.275 



Table 3. Secondary observation detections. Positions are the peak of the emission in the jiggle-maps; note the good agreement with the 
positions predicted by MHW in table|2] Fluxes here are the peak values of the point sources taken directly from the jiggle-maps, with 1-ct 
noise levels estimated according to the SCUBA integration time calculator - the backgrounds are too high in these maps to determine 
noise directly. Fluxes for gel and gc2 have 3-cr upper limits only at 450 fim. For gel and gc2, the spectral index 7 is evaluated for this 
3-<t value (where f oc u" 1 ). For gc6, the spectral index is calculated for the measured values of flux at both wavelengths. 



Source 


R.A. 


Dec. 


850-/xm flux (mjy) 


450-^im flux (mjy) 


7 


gel 


17 46 43.5 


-28 08 21.1 


335 ± 20 


< 1500 


< 2.3 


gc2 


17 44 19.5 


-29 33 33.1 


100 ± 20 


< 1500 


< 4.2 


gc6 


17 43 10.2 


-29 51 47.0 


1069 ± 20 


11318 ± 600 


3.7 



shown in Fig.[S] Each of the spectra shown have had a poly- 
nomial baseline fitted and removed, and have been smoothed 
into bins using a 5-point running average. 



5.3 Discussion of the point source detections in 
continuum and lines 

5.3.1 gel 

Source gel was reasonably strongly detected at 850 ^tm in a 
fairly clear field, but was not detected at 450 /im, although 
there is a hint of a peak in the right place. The position and 
flux of this object in the jiggle-map both agreed well with the 
MHW predictions. The central and average molecular line 
spectra of gel both show the expected signature of a dense 
molecular gas cloud in the Galactic centre (i.e. at ulsr ~ 
Okms -1 ), and so we were able to reject it as an extragalactic 
candidate. 



5.3.2 gc2 

Source gc2 was detected in a field with a visibly high back- 
ground, with an apparent slight extension on the northwest- 
southeast axis. Its position again is extremely similar to that 
predicted, but its peak flux is somewhat below the predic- 
tion by MHW. It was not detected in the 450-^im map. The 
HCO + (3-^2) spectra show a more complicated structure. 
The two positive peaks at about and 20 kms -1 probably 
represent emission from two clouds at slightly different po- 
sitions in the Galaxy, sitting on the same line-of-sight. This 
may be the cause of the slight extension of the source in the 
jiggle-map. The MHW technique is not particularly power- 
ful in separating two closely overlapping point sources, as 
demonstrated here 3 . 



3 However, if two peaks are known to be overlapping, it is possible 
to 'tune' the detection algorithms to separate the two. This will 
be described at more length in K. K. Knudsen et al., in prep. 
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Figure 4. 850-/im jiggle-maps of point-source candidates. Three maps show the only 850-fim detections (top left - object gel; top right 
— gc2; bottom left - gc6); the other (object per3) is included as an example of the non-detections. Units are Jybeam - 1 . The uneven 
edges and holes are caused by 'turning off' noisy bolometers. 



5.3.3 gc6 

The strongest source in the candidate list, gc6 was detected 
strongly in the jiggle-maps at both 850 and 450 /jm. Again 
its position was extremely close to that predicted by MHW, 
and its 850- (im peak flux was 30 per cent lower than pre- 
dicted. It is clearly a compact source. The molecular line 
observations strongly confine this also to be a Galactic ob- 
ject. 

Since both 450- and 850-pim detections were achieved 
for this very bright source, it is worth speculating further 
as to its possible identity. Fitting a standard blackbody to 
the 450/850 colours, a plausible value of emissivity /3<2.0 
is found for temperatures > 80 K. However, fitting a single 
temperature to what is obviously a dense object is probably 
not realistic. 

A 1.4-GHz VLA survey of the area (iLisztJ Il992t 
iGrav et al.lll993Tl suggests that this position is within an 
extended structure which has an integrated flux density of 
several hundred mjy at 1.4 GHz. However, no compact radio 



object was identified within 1.3'of gc6, suggesting an upper 
limit of 3 mjy at 1.4 GHz. There is also an IRAS object 
2'away, but this also seems unlikely to be the true counter- 
part. We are continuing to investigate this unusually bright 
and compact object. 

In summary, the follow-up SCUBA observations demon- 
strate that of 21 objects detected by the MHW routines, only 
3 were genuine point sources. The heterodyne observations 
suggest that all 3 of these are Galactic. The remaining 18 
sources appear to be spurious. 

5.4 Understanding the spurious point sources: 
Altering in scan-maps 

In order to understand why MHW appeared to find so many 
spurious sources, we investigated in detail the original scan- 
map reduction, which revealed an important clue: in the final 
data reduction stage, where the six maps are recombined in 
Fourier space, an option to filter the maps with a low-pass 



8 V.E. Barnard et al. 



00 - 




48 44 17:46:40 36 



Right oscension 




24 17:44:20 16 12 



Right ascension 




14 17:43:10 06 

Right oscension 



Figure 5. 450-^tm jiggle-maps of gel, gc2 and gc6. Units are Jy 
beam - 1 . 
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Figure 6. Spectra of sources gel, gc2 and gc6 (top to bottom). 
All observations were done on 26/08/01. 

filter had been used in both the Galactic centre and Perseus 
maps. The effect of this standard Surf option can be seen 
in Fig. [3 where the power spectra obtained from performing 
the recombination with and without the filtering option for 
the same field are shown. 

The filtering cuts in most steeply at a scale of k > 45, 
which corresponds to a spatial scale of ~ \| , where A is 
the observing wavelength (850 fim) and D is the diameter 
of the dish of the JCMT (15 m). (The 850-^tm beam size 
of 14"corresponds to k ~ 19.) However, the difference plot 
demonstrates that the filtering has an impact at all scales. 
The original intention of the filter was to remove noise which 
appeared on scales smaller than the beam, since this would 
usually be assumed to be caused by instrumental factors. 
However, the information on these small spatial scales is 
important when attempting to detect point sources. In par- 
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Figure 7. Power spectra of filtered and unfiltered maps. The top 
plot shows the two power spectra - the filtered map is represented 
by the solid line and the unfiltered map by the dot-dash line. The 
bottom plot shows the absolute value of difference between the 
two power spectra. 



6 SIMULATIONS 

6.1 Simulating detections of point sources in 
Altered scan-maps 

The scan-map filtering option made a large impact on the 
MHW results by creating spurious sources. However, in the 
earlier stages of analysis, some candidates detected by MHW 
in the filtered maps were ruled out because of their align- 
ment with the positions of known Herbig-Haro objects or 
other sources in the Galaxy. Also, three sources identified 
by MHW were confirmed in the later jiggle-maps. Hence it 
is clear that some real point sources were detectable in the 
filtered maps. Thus we have carried out simulations to deter- 
mine the sensitivity and hence completeness of the analysis 
already carried out, and what selection effects, if any, the 
filtering option had on the detectability of the real point 
source population. 

To evaluate this, simulated point sources were added to 
each of the Galactic centre and Perseus region squares at 
known positions and with known fluxes, and the resulting 
images were fed into the MHW routines following exactly 
the same procedure as used in the original candidate selec- 
tion. The added point sources were spread far enough apart 
that they did not impact upon each other. In a data array 
with the same dimensions as the final map, delta functions 
of unit flux were placed at chosen positions. This array was 
then convolved with a dual Gaussian beam to represent the 
signature of a point source in the uncombined maps. This 
filtering process was repeated six times to produce maps rep- 
resenting each of the combinations of chop throw and direc- 
tion in the uncombined maps. At this stage the six simulated 
source maps were multiplied by a factor to set the flux of the 
point sources, and then each source map was added to the 
respective real uncombined map. These six maps were then 
recombined, both with and without filtering. The fluxes of 
the simulated sources were each gradually increased until 
the MHW detection of each source occurred. 



ticular, noise spikes in a scan-map are smoothed by the fil- 
tering to the right spatial scale to create false positives in 
the MHW algorithms. Whilst of course the initial data re- 
duction included spike removal, the sensitivity of the MHW 
detections went to a deeper level than anticipated at that 
time. 

The maps containing the spurious sources were all re- 
reduced without the filtering option applied, and in every 
case the point source was no longer detected, using the same 
detection parameter requirements. It appears that at each 
location, a noise peak was sufficiently altered by the filtering 
that the MHW routine could not distinguish it from a true 
point source. In the unfiltered maps, each spurious source no 
longer appeared Gaussian to the MHW routine (nor, usually, 
to the eye), whilst the three real sources were still detected. 

These investigations therefore demonstrate that on un- 
filtered SCUBA scan-maps, the MHW routines are a reliable 
detector of point sources. However the filtering option used 
as default in many scan-map images of extended sources will 
mislead detection algorithms into detecting spurious point 
sources. 



6.1.1 Overall sensitivity 

The overall sensitivity of the MHW routines to point sources 
in the filtered maps of the Galactic centre and the Perseus 
star-forming region is shown in Fig. |5] This figure demon- 
strates that the sensitivity in the Perseus region was far 
higher than in the Galactic centre, which is unsurprising con- 
sidering the higher background flux and more complicated 
structure in the Galactic centre maps. In fact the Galactic 
centre count almost seems to have two slopes, with a break 
at around 0.4 Jy. This is probably caused by the variety 
in the background in the Galactic centre, ea sily observed 
by eye (see figures in IPierce-Price et al]l200oTl : those fields 
closer to the Galactic plane suffer from higher backgrounds 
than those in regions slightly offset from the Galactic plane. 
Hence the simulated sources in the squares nearer the Galac- 
tic plane can only be detected at higher fluxes. 



6.1.2 Flux accuracy 

In Fig. [5] the input and output fluxes for simulated point 
sources in the two survey regions are shown. Included on 
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Figure 8. Areal sensitivity of the survey. Area within the Perseus 
region is shown by the dot-dash line, while area within the Galac- 
tic centre region is by the dashed line. The total area is shown 
by the solid line. Effective area here is the fraction of sources de- 
tected with a given input flux, multiplied by the total survey area 
of 4.638 square degrees 



each plot is a line fitted by a least-squares method to eval- 
uate the overall gradient and offset. 

For a line y — mx + c, the values found for the Perseus 
region data were m = 0.80 ± 0.05, c = 0.025 ± 0.005, 
whilst in the Galactic centre data the appropriate values 
were m = 0.85 ± 0.05, c = 0.050 ± 0.005. Hence it ap- 
pears that, on average, the true fluxes are generally slightly 
under-estimated by MHW by 10 or 15 per cent, although 
each graph has a positive intercept. This trend for under- 
estimation may be due to the imperfect beam shapes of 
the scan-map sources (see section 13.21 . The scan-map re- 
construction means that some of the flux of these simulated 
sources was lost into the artifacts surrounding the central 
Gaussian. The MHW technique is not sensitive to these ar- 
tifacts, and so will underestimate the flux of a source. 

The absolute scatter in these graphs appears to stay 
the same at all input fluxes, except perhaps at the very low- 
est fluxes which appear to be systematically over-estimated. 
This suggests that the relative accuracy is worse at lower 
fluxes. This is confirmed in Fig. 1101 where the Relative Flux 
Accuracy (RFA) is shown in comparison with the input 
fluxes required to detect a source. Here the RFA is defined 
as 

Input flux — Detected flux . . 

RrA = . (8) 

Input flux 

Hence values of RFA ~ represent accurate detections, 
negative values represent over-estimations of flux and posi- 
tive values of RFA represent under-estimations. 

The two graphs show similar trends, indicating that the 
fluxes of faint sources tend to be over-estimated, whereas 
brighter fluxes suffer from the effect of the non-Gaussian 
beam as discussed above. At low fluxes this can be under- 
stood as the impact of the noise in the region of the point 
source - a faint source is more likely to be detected if it sits 
in a positive noise region. Whilst the MHW routines will not 
be affected by noise on scales much larger than the beam, 
peaks or holes on scales similar to the beam can affect the 
detection accuracy. 




0.2 0.4 0.6 0.8 1 1.2 



Input flux (Jy) 

Figure 9. Flux accuracy in source detection. The quantisation 
in input flux reflects the resolution of the simulations in terms 
of input flux values tested. The upper plot shows the results for 
the Perseus region fields, while the lower plot shows the Galactic 
centre results. 

6.1.3 Positional accuracy 

The positional accuracy of the detections can also be quanti- 
fied. Fig. llll shows the positional accuracy of each detection 
against the input flux required for first detection of a source. 
It shows broad distributions for both the Perseus and Galac- 
tic centre regions. In these graphs the positional accuracy is 
the absolute distance between the input and output posi- 
tions, in units of 3"pixels. The positional accuracy of the 
detections is extremely reliable with respect to the size of 
the beam. 

6.2 Comparison with unfiltered maps 

To understand the impact of the filtering on the appear- 
ance of true point sources, the same simulations were run 
but without filtering the map at the recombination stage 
in the data reduction. The same data analysis was com- 
pleted, and Fig. 1121 compares the cumulative counts for each 
version of the maps. This figure demonstrates that, in the 
overall population, the sensitivity of the MHW technique is 
not affected by the filtering option. Fig. 1131 shows the de- 
tailed results for the Galactic centre, confirming the overlap 
in detection properties of filtered and unfiltered maps. It is 
clear that the filtering procedure may create spurious extra 
sources but does not lose real ones. 
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Figure 10. Relative flux accuracy in source detection. Lines are 
fitted to distribution — see section 16.1.21 Upper and lower plots 
are as in Fig. I§1 



Figure 11. Positional accuracy in source detection. The quan- 
tisation here reflects the pixelised nature of the maps, where 1 
pixel = 3". Upper and lower plots are as in Fig. [5] 



6.3 Possible problems with MHW and SCUBA 
scan-maps: Loss of point sources 

During the analysis of the simulations, some strange be- 
haviour in the detection history of individual sources became 
apparent: a simulated source detectable at a lower input flux 
was sometimes not detected at the next higher input flux. 
Also, real sources in the fields sometimes disappeared from 
the detection lists as the fluxes of the additional simulated 
sources increased. Not every source nor every field was af- 
fected however. It appears that two modes of 'point source 
loss' were occurring. 

6.3.1 Mode A: Increasing \ 2 as fl ux increases 

In this mode, point sources were not detected by the MHW 
routines because their x 2 values were too high. The \ 2 P a ~ 
rameter fsection l4.2H is a measure of how similar the shape of 
a point source is to a Gaussian. The limiting value of 4 used 
here was determined from previous simulations of Planck 
detections. Whilst the Planck simulations were run for the 
same wavelength as these SCUBA observations, the simu- 
lated Galactic emission was of a general cirrus component 
on larger spatial scales, and so was not representative of the 
high-flux and variable structures in these maps. The Planck 
beam was also assumed to be perfectly Gaussian. The beam 
map in Fig. [5] shows that the beam of the SCUBA scan- 
map mode is not a perfect Gaussian; in particular there are 
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Figure 12. Unfiltered simulation results: cumulative counts. The 
overall results for the filtered sources are shown by the solid line 
with crosses and the same results for unfiltered sources are shown 
by the dashed line with solid circles. 



features related to the image reconstruction technique. As 
the flux of a simulated source increases, imperfections in 
the beam become more important relative to the local back- 
ground, increasing the \ 2 value of the source and preventing 
its detection. 

To check this explanation, the simulations were re- 
peated with two variations. Firstly, a single point source 
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Figure 13. Comparison of filtered and unfiltered results for the 
Galactic centre. Filtered sources are shown as squares, as ear- 
lier. Unfiltered sources are shown as solid circles. Results for 
the Perseus region showed similar overlap but are not shown for 
brevity. 



was added to a map, rather than several at once, to check 
that the beam profiles of adjacent simulated sources were 
not interfering. The detection history for a particular source 
was found to be identical whether it was added alone to a 
map or multiple sources were added at once, indicating that 
overlapping of sources was not an important factor in these 
simulations. Secondly, the individual sources were replaced 
by perfect Gaussians of the same height and beamwidth. In 
this case no loss of sources was observed at all. These two 
results confirm that it is likely to be the imperfectly Gaus- 



sian beams of the individual scan-map point sources which 
causes their loss at higher fluxes. 

Having discovered this effect in the simulations, it now 
seems likely that this is the cause of the very weak detection 
of the very bright source gc6 in the initial maps (see section 
14.21 . This exceptionally bright source, which from the jiggle- 
map observations appears still to be a real point source, was 
only detected by the MHW routines with some relaxation of 
the detection parameters, to allow a higher \ 2 value of 5.26. 
Given its flux, its seems probable that its high \ 2 value was 
due to the impact of the outer imperfections in the beam. 

In general this mode of point source loss affects only the 
brighter sources in SCUBA scan-maps. A check by eye of all 
fields being surveyed identified any possible problem sources, 
and we are confident that no real bright point sources were 
lost in this manner, gc6 having been found. The flux level 
at which sources are lost varied from position to position, 
though generally this loss mode only occurred at input fluxes 
> 800 mjy; it depended on the local features of the back- 
ground around the source, since these determined when the 
imperfections of the beam become significant. 



6.3.2 Mode B: the contribution of point sources to 



.(«) 



A second mode of point source loss was observed in the sim- 
ulation results. Unlike the previous mode, this could occur 
at both low and high flux levels, and affected the real point 
sources in the map, as well as the simulated sources. 

In this loss mode, sources fail the second detection 
requirement, that the signal-to-noise in wavelet space, 
D w (Ro V %), is > 5 which was also set from the Planck sim- 
ulations. An implicit assumption within this formalism is 
that the point sources themselves are not affecting the value 
of o~ Wn (Ro P t), the noise in wavelet space. a Wn (R op t) is cal- 
culated from the properties of the whole map, and ideally 
should be determined from a background without any point 
sources. In practice this is not plausible, as the Galactic cen- 
tre and Perseus backgrounds are hard to model. Instead, the 
point sources in a map can contribute to the value of a Wn (R) 
at high enough flux levels, and so affect their own chance of 
detection. 

Similar checks were carried out to test the behaviour 
of this mode of point source loss as for mode A. The effect 
was reduced but did not completely disappear with perfectly 
Gaussian point sources. Also the number of point sources 
added to a map did make a difference to the results. This 
mode of point source loss is related therefore to the overall 
properties of the image map used, as well as to the properties 
of the sources. In smaller maps the presence of bright point 
sources is more likely to make an impact on the value of 
a Wn (R op t), possibly to the extent that fainter point sources 
will not be detectable. For regions where the background 
cannot be simulated, two ways to avoid this effect would be 
to either use larger maps, or to alter the detection require- 
ments to lower D m (R op t) values for maps which are either 
crowded or contain bright sources, in order to be sure of de- 
tecting the fainter sources. The latter of course will increase 
the number of spurious detections as well though. This mode 
of point source loss will only have affected the survey results 
for rare fields in which there were very bright point sources. 
Since each field was checked by eye, we are again confident 
that no real sources were lost in this way. 



The very bright SCUBA galaxy count 13 



6.4 Summary of simulation results 

In general the simulations suggest that the survey carried 
out with the MHW algorithms has produced reliable results, 
and the small average errors in flux and position estimation 
are easily characterisable. It appears that the filtering option 
used in the construction of the final images did not affect 
the detection of real point sources as a population. Two 
mechanisms which may lead to non-detections of real point 
sources have been identified, but it seems unlikely that either 
strongly affected the result of this survey. 



7 FINAL RESULTS 

The final results of this paper are constraints put on the 
counts of SCUBA galaxies at bright fluxes. This is some- 
what simplified (!) by the fact that no extragalactic candi- 
dates were detected. However the results of the simulations 
indicate the sensitivity of the source detection method used, 
and so were used to set appropriate upper limits on the pos- 
sible population. 

Poiss on confidence li mits for small numbers were pub- 
lished by lOehrela (jl98f)fl . The upper limit corresponding 
to a l-<7 Gaussian limit is 1.841. 4 For flux values in the 
range suggested by Fig. |H| the upper limit was calculated 
as 1.841 / ' Atotai- This resulted in the final error bars listed 
in Table [I] Below flux limits of approximately 200 mjy the 
Galactic centre survey is relatively insensitive, whereas the 
Perseus survey is complete above about 300 mjy. Hence the 
l-cr upper limit on the 850-/im count is 53 per square degrees 
at 100 mjy and 2.9 per square degree at 50 mjy. 

In Fig. 1141 the upper limits in table 0] are plotted, with 
two lines. The dashed line is a Euclidean slope (oc S ,_3//2 ), 
with the highest intercept consistent with all the upper limit 
error bars: at S — 130 mjy, the count of galaxies is 1 per 
square degree. The strongest constraints on this normalisa- 
tion are the error bars at 150 - 200 mjy, where the combined 
bright survey becomes sensitive over more than 50 per cent 
of the available area (see Table 2J. The solid line represents 
the power law fitted in Fig. 0to the fainter SCUBA galaxy 
counts. Also plotted are th e results of the two brightest 
blank -field SCUBA surveys JScott et al.ll2002t iBorvs et all 
120031 . see table 1 for details). 

The count of galaxies at fluxes lower than 10 mjy is 
clearly super-Euclidean. However both the 8 mjy survey and 
HDF super-map results appear to show a slight turnover at 
their brightest fluxes, ~ 10 - 15 mjy. Whilst these results 
are of course subject to larger uncertainties than at lower 
fluxes, the counts at these bright fluxes are more consistent 
with the Euclidean slope fitted to the upper limits derived 
in this paper. This supports the idea that the steep turnover 
seen in the ~ 15 mjy counts is real and continues to brighter 
fluxes. 

This r esult also may be consistent with a suggestion by 
lLawrencel i200lf) that a minority of the blank-field SCUBA 
sources, even at fluxes as low as 2 mjy, might be local, very 
cold, dark clouds at distances of around 100 pc. Such sources 
would be expected to show a Euclidean trend at fluxes 

4 The equivalent limits for 2- and 3-cr are 2.303 and 2.996. For a 
zero count there is no lower limit. 
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Figure 14. Euclidean and modelling predictions for counts. The 
dashed line shows the highest possible Euclidean slope allowed 
by these data. The solid line represents the slope of the count de- 
termined at fainter fluxes, as shown in Fig. 1. The large symbols 
at higher fluxes show predictions from si x galaxy evolution mod- 
els predictin g bright counts at 850 fan'. |F^anc^schin^^alJ^993 
(diamonds) lOuiderdoni et ^lj|l99 ^ (circle) ; [ Granato et al .11200 ll 
(crossed circle) ; |D cvriendt &: GuiderdonilpOOO j (curved squares) ; 
feowan-^obingonl l200ll (diagonal crosses) : IPearsonll200ll (square); 
iBlain et al.l 1999 (open crosses). The smaller symbols with error 
bars to the left show the results from the two brightest blank-field 
SCUBA surveys to date (triangles: Scott et al. 2002; stars: Borys 
et al. 2003) 

brighter than about 10 mjy, although recent confirmation 
of cosmological redshifts for the majority of radio-sel ected 
SCUBA-detected objects jChapma.n et alJl200.1 l2004li con- 
firm that any such Galactic population is not numerically 
significant. The Euclidean slope plotted here corresponds to 
about 20 per cent of the count found bv lScott et alJ l)2002l) 
at 8mJy, and nearly 90 per cent at 12 mjy. 

Also plotted on Fig. ll4l are the predictions for the bright 
850-pim count from various galaxy evolution models which 
give bright count predictions. These are all consistent with 
the limits found from the Galactic centre and Perseus sur- 
veys. 



8 CONCLUSIONS 

In this paper an application of the Mexican Hat Wavelet 
to detect SCUBA galaxies was presented. Bright candidates 
were detected in scan-maps of the Galactic centre and the 
Perseus star-forming region. Follow-up continuum and spec- 
tral line observations suggested that none of the 21 initial 
candidates were extragalactic; furthermore, only three ap- 
peared to be genuine point sources. Re-reducing the scan- 
map data without filtering its Fourier spectrum removed the 
spurious point source candidates. 

Simulations of the MHW on the Galactic centre and 
Perseus fields demonstrated that the technique would have 
detected point sources flux over the flux range 50 - 1000 
mjy. The flux and positional accuracies of recovered sources 
were shown to be high. A problem involving the loss of point 
sources was revealed, but overall the MHW technique has 
proved robust, particularly in the presence of very compli- 
cated backgrounds. The two key advantages of the MHW 
techniques are i) the quantitative evaluation of both the 
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Table 4. Final results. Columns 2, 3 and 4 show the area over which the survey was sensitive to the flux limit in, respectively, the 
Perseus region, the Galactic centre and the two regions combined. Column 5 gives the 1 a upper limit as calculated as 1.841/A tota i. 



Flux limit (Jy) 


A pors (sq. deg.) 


A gc (sq. dcg.) 


Atotai (sq- deg.) 


a A (per sq. deg.) 


0.05 


0.028 


0.007 


0.035 


53 


0.07 


0.070 


0.007 


0.077 


24 


0.1 


0.616 


0.018 


0.634 


2.9 


0.15 


1.930 


0.107 


2.037 


0.9 


0.2 


2.657 


0.280 


2.937 


0.6 


0.3 


2.797 


0.736 


3.533 


0.52 


0.4 


2.797 


1.123 


3.920 


0.47 


0.5 


2.797 


1.333 


4.130 


0.45 


0.6 


2.797 


1.440 


4.237 


0.43 


0.7 


2.797 


1.528 


4.325 


0.43 


0.8 


2.797 


1.602 


4.399 


0.42 


0.9 


2.797 


1.727 


4.524 


0.41 


1.0 


2.797 


1.790 


4.587 


0.40 


1.2 


2.797 


1.841 


4.638 


0.40 



shape and flux of the source, and ii) that no assumptions 
are required about the characteristics of the noise in an im- 
age, which is often an implicit feature of other routines. 

For maps made with SCUBA, the MHW technique 
is particularly appropriate when the beam is complicated, 
since the central Gaussian feature to which the MHW rou- 
tines are sensitive will still remain. The MHW is also likely 
to be the most robust technique in maps with high or vary- 
ing backgrounds. This work has demonstrated that point 
sources can be reliably located and measured in SCUBA 
scan-maps when a suitable detection algorithm is used, and 
that the scan-mapping mode should now be considered a 
possible tool for extragalactic submillimetre surveys. 

Using the results of simulations, upper limits have been 
established for the count of SCUBA galaxies at fluxes > 50 
mjy. The Euclidean slope then consistent with these upper 
limits approaches the previous brightest counts at fluxes ~ 
15 mjy, suggesting that the turnover seen in the brightest 
blank-field surveys is real and continues to very bright fluxes. 
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